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Abstract 

We perform a Dalitz plot analysis of the decay Df — > ivr + i^ _ 7r + with the CLEO-c data set of 
586 pb _1 of e + e~ collisions accumulated at yfs = 4.17 GeV. This corresponds to about 0.57 million 
DfD* T pairs from which we select 14400 candidates with a background of roughly 15%. In contrast 
to previous measurements we find good agreement with our data only by including an additional 
/o(1370)7r + contribution. We measure the magnitude, phase, and fit fraction of K*(892)° K + , 
0(1O2O)tt+, i^ *(1430)i^+, / (980)vr + , / (1710)vr + , and /o(1370)tt + contributions and limit the 
possible contributions of other KK and Kir resonances that could appear in this decay. 



* Deceased 
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I. INTRODUCTION 



The decay Df — > K + K~tt + is among the largest known branching fractions for the D s 
meson. For some time the mode — > 0(lO2O)7r + was used as the normalizing mode for D s 
decay branching fractions, typically done by choosing events with the K + K~ invariant mass 
near the narrow 0(1020) peak. Observation of a large contribution from — > / (980)7r + 
[l| makes the selection of — > 0(lO2O)7r + dependent on the range of K + K~ invariant 
mass chosen; the observed yield of non-0 contributions can be larger than 10% 0]. This 
is an unacceptably large uncertainty for a normalizing mode and we proposed [2J that the 
branching fraction for Df — > K + K~tt + in the neighborhood of the peak, without any 
attempt to identify the 07r + component as such, could be used for D s normalization. Relating 
the .D+ — > K + K~n + branching fraction in |2j to the rates for such phase space-restricted 
subsets requires an understanding of the resonance contributions to the final state. The only 
published Dalitz plot analysis [3j has been done by E687 pj using 701 signal events. The 
FOCUS Collaboration has studied this decay in a Dalitz plot analysis in an unpublished 
thesis 0] and a conference presentation j^]. 

Here we describe a Dalitz plot analysis of Df — > K + K~n + using the CLEO-c data 
set which yields a sample of over 12, 000 signal candidates. Charge conjugation is implied 
throughout except where explicitly mentioned. The next section describes our experimental 
techniques, the third section gives our Dalitz plot analysis formalism, the fourth describes 
our fits to the data, and there is a brief conclusion. 

II. EXPERIMENTAL TECHNIQUE 

CLEO-c is a general purpose detector which includes a tracking system for measuring 
momenta and specific ionization of charged particles, a Ring Imaging Cherenkov detector to 
aid particle identification, and a Csl calorimeter for detection of electromagnetic showers. 
These components are immersed in a magnetic field of 1 T, provided by a superconducting 
solenoid, and surrounded by a muon detector. The CLEO-c detector is described in detail 
elsewhere jsj. 

We reconstruct the Df — > K + K~7t + decay using three tracks measured in the tracking 
system. Charged tracks satisfy standard goodness of fit quality requirements Q. Pion and 
kaon candidates are required to have specific ionization, dE/dx, in the main drift chamber 
within four standard deviations of the expected value at the measured momentum. 

We use two kinematic variables to select — > K + K~n + decays, the candidate invariant 
mass 

m inv = m(K + K~7v + ) or Am inv = m inv - m Ds , (1) 
and the beam constrained mass 

m BC = \f ^bcam - Pd oi Att^c = m BC - m BC {D s ), (2) 

where 771^=1968.2 MeV/c 2 [8j is the D s mass, -Ebeam is the beam energy, p D is the mo- 
mentum of reconstructed candidate, and m-Qc(D s ) = 2040.25 MeV/c 2 is the expected 
777bc value of the D s meson in the process e + e~ — > D*D S at = 4170 MeV. We require 
|Am inv | < 2cr(m inv ), |AmBc| < 2<r(mBc), where the resolutions cr(m inv ) = 4.8 MeV/c 2 
(4.79±0.05 MeV/c 2 in single Gaussian fit), and a(m BC ) = 2 MeV/c 2 (1.89±0.02 MeV/c 2 ) 
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represent the widths of the signal peak in the two dimensional distribution. When there are 
multiple D s -meson candidates in a single event we select the one with smallest Amec value. 

We use a kinematic fit to all 3-track combinations which enforces a common vertex and 
mass constraint. The kinematic fit-corrected 4- momenta of all 3 particles are used to 
calculate invariant masses for further Dalitz plot analysis. The resolution on the resonance 
invariant mass is almost always better than 5 MeV/c 2 . 

After all requirements, we select 14400 candidate events for the Dalitz plot analysis. 
The fraction of background, 15.1%, in this sample is estimated from the fits to the m- mv 
spectrum shown in Figure [H In most fits to the Dalitz plot we constrain the value of the 
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FIG. 1: The 

m inv distribution. The vertical (online blue) lines show the ±2<r signal region. 



signal fraction, / sig = 84.90±0.15%. In cross checks we use a set of sub-samples, splitting the 
data by time of observation and by sign of Z) s -meson charge, Df and D~ . We also consider 
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samples with tight (lxl standard deviations in mec and m; nv ) and loose (3x3 three 
standard deviations) selection versus the standard selection, as well as samples of D s mesons 
produced in D* — > D s ^ decays, selected with a displaced signal box using m B c low band 
(|m BC -2025 MeV/c 2 | < 4cx(m BC )) and m BC high band (|m BC -2060 MeV/c 2 | < 4cr(m BC )). 

To determine the efficiency we use a signal Monte-Carlo (MC) [9] simulation where one 
of the charged D s mesons decays in the KKtt mode uniformly in the phase space, while 
the other D s meson decays in all known modes with relevant branching fractions. In total 
we generated 10 6 Df and D~ signal decays. These underlying events are input to the 
CLEO-c detector simulation and processed with the CLEO-c reconstruction package. The 
MC-generated events are required to pass the same selection requirements as data selected 
in the signal box. We only select the signal-side D s mesons which decay uniformly in the 
phase space, separating them by charge. 

We analyze events on the Dalitz plot by choosing x = m 2 (K + K~~) and y = m 2 (K~7T + ) 
as the independent (x,y) variables. The third variable z = m 2 (K + n + ) is dependent on x 
and y through energy and momentum conservation. We do not expect any resonant sub- 
structure in the K + 7r + invariant mass; with these Dalitz plot variables any structure in z is 
due to reflections of structures in x and y. Figure [2] shows the Dalitz plot. Besides the clear 
0(1020) and K*(892) signal, no other narrow features are clearly observed. The variation of 
the population density along the resonance band clearly indicates that these resonances are 
spin one as the amplitude for a spin-one resonance should have a node in the middle of its 
band. There is a significant population density in the node region of the 0(1020) resonance, 
indicating that there is likely to be an additional contribution. 

To parametrize the efficiency, e(x,y), we use a third-order polynomial function with 
respect to the arbitrary point (x c , y c ) = (2,l) (GeV/c 2 ) 2 on the Dalitz plot times threshold 
functions in each of the Dalitz variables to account for the loss of efficiency at the edges of 
the Dalitz plot, such that 

e(x,y) = y)T(x)T(y)T(z(x,y)). (3) 

With x = x — x c and y = y — y c , the efficiency is the product of the polynomial function, 

£ P oiy(x,y) = l + E x x + E y y + E x 2X 2 + E y2 y 2 + E x sx 3 + E yi y 3 + E xy xy + E x 2 y x 2 y + E xy 2xy 2 , (4) 

For each Dalitz plot variable, v (= x, y or z) the threshold function is sine-like with 

T ( v \ = i l E c,v + (1 - E c>v )] x sm(E th , x \v - v max \), at < E th>v x \v - w max | < tt/2, 
' \ 1, at E thjV x\v- v max | > vr/2, 

(5) 

All polynomial coefficients, E x , E y , E x 2, E y2 , E x 3, E y a, E xy , E x 2 y , E xy 2, E CjV , and E t h, v are 
fit parameters. Each variable v has two thresholds, t> m in and f max - We expect low efficiency 
in the regions v ~ v max only, where one of three particles is produced with zero momentum 
in the D s meson rest frame and thus has a small momentum in the laboratory frame. 

The simulated signal sample is used to determine the efficiency. Table [I] shows the results 
of the fit to the entire signal MC sample of — > K + K~tt + events selected on the Dalitz 
plot. The polynomial function with threshold factors describes the efficiency shape very 
well for our sample. We also fit separately the signal MC sub-samples for — > K + K~tt + 
and D~ — > K~K + tc~ decays, for simulations of early and late datasets, and for tight and 
loose signal boxes. In cross-checks with sub-samples we fix the threshold parameters to their 
values from the central fit in order to remove correlations with other polynomial parameters. 
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FIG. 2: The Dalitz plot for the data. 

We find that the variation of the efficiency polynomial parameters is small compared to their 
statistical uncertainties. In fits to data we use this efficiency shape with fixed parameters, 
and constrained variation is allowed as a systematic check. 

The shape for the background on the Dalitz plot is estimated using data events from a 
777.bc sideband region, \m-Qc — 1900MeV/c 2 | < 5cr(mBc)- We only consider events from the 
low mass m BC sideband as the high mass sideband is contaminated by signal events due to 
initial state radiation. To parametrize the background shape on the Dalitz plot we employ 
a function similar to that used for the efficiency, shown in Eq. HJ We add incoherently to 
the polynomial two peaking contributions to represent If* (892) and 0(1020) contributions 
described with Breit-Wigner functions with floating normalization coefficients, Bk* and B^,, 
respectively. Figure [3] and Table [TT] show results of the fit to the background polynomial 
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TABLE I: Fit parameters for describing the efficiency across the Dalitz plot. 



Parameter 


Value 


TP 


n noo i r\ mo 

U.U2o±L).Ulz 


TP 

£jy 


U.Uo 1 ±U.U14 


TP 


— U.oU I ±U.U14 


TP 


— U.0ZD±U.Uo4 


TP 

Ey2 


n oni _i_n hq/I 


TT 

&xS 


U.zoz±U.Uzb 


TP 


U.953±U.U78 


TP 

&xy2 


U.887±U.U98 


Ey 3 


0.004±0.05f 


Eth,x 


3.23±0.18 


Eth,y 


2.53±0.13 


Eth,z 


2.61±0.13 


E 

J -'c,x 


0.166±0.042 


E 

&c,y 


0.320±0.034 


E c , z 


0.338±0.034 



TABLE II: Fit parameters for the background sample. Values in parentheses show an uncertainty 
or variation of the last significant digits. 



Parameter 


Value 


B x 


-0.23±0.11 


By 


0.06±0.13 


B X 2 


-0.29±0.12 


B X y 


-0.99±0.29 


By2 


-0.47±0.32 


B X 3 


0.77±0.23 


B x 2y 


1.98±0.67 


B X y2 


2.24±0.84 


By3 


0.56±0.47 


Bff, 


0.000f6f(23) 


Br* 


0.00f44(28) 



function for our sample. We also consider the variation of the background shape parameters 
for sub-samples, split for Df and Dj, for earlier and later datasets, and for tight and loose 
cuts on background selection box. The variation of the shape parameters is small compared 
to their statistical uncertainties. Furthermore, in fits to data we use the background shape 
with fixed parameters, and constrained variation is allowed as a systematic cross check. We 
also allow the size of the narrow resonance contributions to the background to float freely 
as a systematic variation. 
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FIG. 3: Projections of the fit to the background shape described in the text, line, displayed over 
the data, dots, in the background box. 



III. FORMALISM 



This Dalitz plot analysis employs the techniques and formalism described in Ref. |10j that 
have been applied in many other CLEO analyses. We use an unbinned maximum likelihood 
fit that minimizes the sum over N events: 



A' 



£ = -2^1ogP(x n ,y n ), (6) 



n=l 



where V(x, y) is the probability density function (p.d.f.), which depends on the event sample 
being fit, 

Af e e(x,y) for efficiency; 

V(x,y) = \ AfBB(x,y) for background; (7) 

fsi g Ns\M{x,y)\ 2 e{x,y) + (1 - f sig )N B B(x,y) for signal. 
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The shapes for the efficiency, e(x,y), and background, B(x,y), are discussed in the pre- 
vious section. The signal p.d.f. is proportional to the efficiency-corrected matrix element 
squared, \Ai(x,y)\ 2 . As described above, the signal fraction, / sig , is defined from the in- 
variant mass spectrum. The background term has a relative (1 — / S i g ) fraction. The effi- 
ciency, signal, and background fractions are normalized separately, 1/Af £ = J e(x,y)dxdy, 
1/A/s — / \-M.(x, y)\ 2 e(x, y)dxdy, I/Mb = / B(x,y)dxdy, which provides the overall p.d.f. 
normalization, / V(x,y)dxdy = 1. The matrix element is a sum of partial amplitudes, 

M = ^c R xW R xn R x^x Ti, (8) 

R 

where Wr depends on the spin of resonance R. The factor Qr is the angular distribution 
for the resonance, and the factors J 7 ^ and are the Blatt-Weisskopf angular momentum 



barrier-penetration factors In our standard fit the complex factor cr = aRt l ^ R is 

represented by two real numbers, an amplitude clr and a phase <pR. These are included in 
the list of fit parameters and can be left to float freely or fixed. 

Assuming the decay chain d — > Rc — > abc we may write the angular distribution, 

nr = i. 

«r = k - < + (mg - m;2)( r° <9) 



m ab 



2 

ab ' m ab 

where is the mass of the decaying particle and m a , mb and m c are the masses of the 
daughters; m a b, m ac , and m^c are the relevant invariant masses. These expressions for angular 
distributions can be obtained from covariant-tensor formalism or from orbital momentum 
partial waves decomposition using Legendre polynomials Pl (cos 9), where 9 is the angle 
between particles a and c in the resonance R rest frame. 

For regular resonances such as IT" (892), 0(1020), IT* (1410), fQ(1430), etc., we use the 
standard Breit-Wigner function, 

W R {m) = — ^ (10) 

multiplied by the angular distribution, Ql, and the Blatt-Weisskopf form factors T^{q) and 
F%(q) for the D -meson and resonance R decay vertexes, respectively. We assume that the 
mass dependent width has the usual form 

T{m) = r*^(^) 2i+ Vl(P x r R )}\ (11) 

where P is the decay products' momentum value in the decaying particle rest frame and tr 
is the effective resonance radius. The form factors T^q) and T^q) in Eqs. [81 and [TT1 are 
defined in the Blatt-Weisskopf form [Hj 

L = 0: J$(q) = l, (12) 
L = l: ^(g)=A^x [l + g 2 ] _V2 , (13) 

-1/2 



L = 2: ^(g)=A^x 9 + 3 9 2 + 9 4 , (14) 
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where the label V stands for D or R decay vertex, q = P x ry, ry is an effective meson 
radius, and My is a normalization constant defined by the condition jFy(P^ x ry) = 1, 
where Pr is the products' momentum value at m = tur. 

The Wr parameterization of the / (980), whose mass, mj , is close to the KK production 
threshold, uses the Flatte formula 



WflH = , ^ 7^ , : j-z (15) 



///?, - m- - i J^ab 9RabPab(m) 



where ab stands for n°ir , n + n , K + K , and K°K°, and p a b( m ) = 2P a /m is a phase space 
factor, calculated for the decay products momentum, P a , in the resonance rest frame. We 



use the following isospin relations for the coupling constants <7/ 07r +7r- = w2/3g/ 07r7r , gf 07T o n o 



1/3^/0^, and 9 hK ^ = 9f K+K- = yU 2 9f KK- Tneir values, shown in Table Hill are 
taken from the BES experiment 13 



We model a low mass K + 7r S wave, also known as k or If (800), using a complex pole 
amplitude proposed in Ref. [141 ]. 



w«M = ( 16 ) 



m 2 . — 



where m K = (0.71 — z0.32) GeV is a pole position in the complex s = m 2 (K + 'K~) plane 
estimated from the results of several experiments. 

In this analysis we use or test all known K~tc + and K + K~ resonances recognized by 
the Particle Data Group (PDG) [§| which can be observed in the phase space of the — > 
K~K + -k + decay. These are listed in Table IIHI One could expect a contribution in the 
K + K~ mass spectrum from the /o(980) and ao(980) scalar resonances. Their K + K~ mass 
spectra have similar, but not well defined shapes. If both amplitudes are allowed to float 
simultaneously in the fit, they show a huge destructive interference, which is sensitive to 
their shape parameters. The /o(980) contribution dominates jl] in the — > 7r + 7r + 7r~ decay, 
which has a large branching fraction, B(Df — > 7r + 7r + 7r~) = (1.22 ± 0.23)%. The relevant 
coupled channel of the ao(980) has not been observed in the — > i]7i°7i + decay. In this 
analysis we consider the /o(980) contribution only. 



IV. FITS TO DATA 

First, we analyze our data with the model used by E687 [l[. Their isobar model contains 
five contributions, K*(892)°K+, 0(1O2O)tt+, ^(1430)^+, /o(980)tt + , and / (1710)vr + . In 
our analysis of D + — > K~7t + 7t + and D + — > K~K + it + decays we find a if* (892) width that 
is smaller than the world average value from the PDG 0. Thus we let the mass and width 
of if* (892) float in the fit. Results are shown in Table HVl In this table and all succeeding 
tables, the units of the amplitudes are arbitrary (a.u.). We find that the sign of the 0(1020) 
contribution is opposite to the sign obtained by E687, but all other results are consistent 
within quoted uncertainties. We find that this fit to our data sample has a poor x 2 / u i 
where v is the number of degrees of freedom, giving a very small fit probability. The x 2 
is calculated over adaptive bins, similar to our previous analysis [3]. This model does not 
represent our data well especially in the range of 1.1 < m 2 KK < 1.5 GeV 2 /c 4 . 

The E687 model contains five resonances. Two of them, K * (892) and 0(1020), are clearly 
seen on the Dalitz plot. The other three, ^(1430), /o(980), and /o(1710), are too wide to 
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TABLE III: Parameters of contributing resonances. 



Resonance J PC 


Mass (MeV/c 2 ) 


Width (MeV/c 2 ) 


Kit states 


if* (892) 


l - 


896.00±0.25 


50.3±0.6 


K*(UW) 


1" 


1414±15 


232±21 


Kg; (1430) 


0+ 


1414±6 


290±21 


iff (1430) 


2+ 


1432.4±1.3 


109±5 


if* (1680) 


1- 


1717±27 


322±110 


K 


0+ 


SRm=710 


9 m=-310 


K+K~ states 


/ (980) 


0++ 


965±10 


g n7T =406 
g K K=800 


a (980) 


0++ 


999±1 


g V7T =620 
g K K=500 


0(1020) 


1— 


1019.460±0.019 


4.26±0.05 


/ 2 (1270) 


2++ 


1275.4±1.1 


185.2l|i 


a 2 (1320) 


2++ 


1318.3±0.6 


107±5 


/o(1370) 


0++ 


1200 to 1500 


200 to 500 


a (1450) 


0++ 


1474±19 


265±13 


/o(1500) 


0++ 


1507±5 


109±7 


/a (1525) 


2++ 


1525±5 


73^ 


/o(1710) 


0++ 


1718±6 


137±8 


0(1680) 


1~ 


1680±20 


150±50 



be easily discerned. To check their significance we remove them one-by-one from the total 
amplitude and check the fit results. In all fits where we remove one resonance the fit quality 
is degraded, increasing \ 2 / v by more than 0.6, compared to our central fit. Thus, we assume 
that all five resonances from E687 model are significant. 

In order to get better consistency between the model and data, we try to improve the 
E687 model by adding contributions from the other known resonances listed in Table IIHI 
The results of these fits are shown in Tables |V] and IVII as a variation of the fit parameters 
with respect to the central case. In all cases the fit quality is improved and each additional 
resonance has a significant magnitude. We conclude that the five resonance model based on 
E687 results does not fully describe the data sample. The largest fit quality improvement 
is achieved in the case of additional S-wave contributions: /o(1370); non-resonant (NR); 
a (1450); and k. Adding the / (1370) contribution gives the largest improvement of the fit 
quality, A% 2 = — 100 for two fewer degrees of freedom. 

We consider a six-resonance model, called Model A, containing K *(892)°if + , 0(1O2O)7T + , 
K*(U30)K+, /o(980)7r+, / (1710)vr+, and /o(1370)tt + contributions. Model A is simply the 
E687 isobar model with an additional /o(1370)7r + contribution. Results with this model and 
fit projections are shown in Fig. HI We repeat the previous procedure and include one-by- 
one additional resonance and check the significance of its parameters and consistency of the 
p.d.f. with our data sample. Results are shown in Tables IVII I and IVIIII For Model A we 
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TABLE IV: Comparison of CLEO-c results with E687 using the E687 isobar model. Shown are the 
fitted magnitudes, a in arbitrary units, the phases (<fi) in degrees, defined relative to the K*(892)°7r + 
amplitude, and the fit fractions (FF). 



Mode 


Parameter 


E687 


CLEO-c [PDG1 




Q 


^iixeti j 


i ^11-x.eu. j 




6 (°) 


n (fixed") 

\J I 11 jVv3 v_l / 


ffixed"l 




m (MeV/c 2 ) 




895 8+0 5 [896 00±0 251 




r fMpV/r 2> i 

L I Ivlc V / C J 








FF f%1 


47 8±4 6+4 


48.2±1.2 




(I 


N/A 


1 76+0 12 




(h (°\ 
9 \ I 




1 4.CL_|_C 

± _l_ O 




FF (%) 


9.3+3.2+3.2 


5.3+0.7 


c6fl020W+ 





N/A 


1 1 5+0 02 




9 \ I 


1 7o_|_on_i_o z i 

± 1 ill ZjI 1 1 Z.'-f 


-1 5+4 




FF (%) 


39 6±3 3+4 7 


42.7+1.3 


fn(980)7r + 


a 


N/A 


3.67+0.13 




*n 


159+22+16 


156+3 




FF (%) 


11.0±3.5±2.6 


16.8+1.1 


/o(1710)7T+ 


a 


N/A 


1.27+0.07 




<M°) 


110+20+17 


102+4 




FF (%) 


3.4+2.3+3.5 


4.4+0.4 


£ FF (%) 




111.1 


117.3+2.2 


Number of events on DP 






14400 


Number of Signal events 




701±36 


12226+22 


Goodness 


X 2 /" 


50.2/33 


278/119 



do not find any additional resonances with significant magnitude, the fit quality does not 
significantly improve, and thus we take this model for our central result. For each additional 
resonance we estimate an upper limit on its fit fraction at 90% confidence level, as also shown 
in Tables IVIII and IVIIII We conclude that the six-resonance Model A p.d.f. gives a good 
description of our data sample. 

For Model A we test the resonance shape parameters by floating the mass and width, or 
two coupling constants in case of /o(980), for each resonance. Results of these fits are shown 
in Tables IIXI and |2D We find that all parameters are consistent with their central fit values 
used in the fit with Model A. 

To estimate systematic uncertainties of the fit parameters, we apply numerous variations 
to the fitting procedure and look at the change of the fit parameters from the central result. 
We consider sub-samples where the data is split into earlier and later datasets, and D~ 
decays, and selected using tight and loose signal boxes. These are shown in Table IXII These 
results are obtained with fixed parameters for efficiency and background functions from 
Tables H] and [TTl We also consider fits with floating efficiency or background parameters in 
Table IXIII In these fits all polynomial coefficients for the efficiency or background including 
resonance background amplitudes float freely, but we fit simultaneously two samples of 
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TABLE V: Fits to CLEO-c data using the E687 model with additional K ir + resonances. For the 
contributions that do not change the entries in the table are changes from the E687 model. 



Pqvq t~v~i py 
± cLL dint; LCI 




NR 


K*(1410) ET 2 *(1430) iT"(1680) 


t- 

rt 


m K*(892) 


one; o_i_n K 


n n 


-n zl 


— n i 

U. 1 


— 1 9 


n q 


1 if* (892) 


AA 9_i_i n 




-1.3 


0.3 


1 
— Z. 1 


-D 3 


fl r , > _ . „ x foil \ 

"^(1430) [O..U.) 


i 7fi-L.n 1 9 

1 . i UZIZU. 1^ 


1 . IU 


-0.02 


0.14 


o 


U.(JO 


9 'K* (1430) I ) 




-4 9 

■4:. 


4 


7.3 


-4 


7 


u /o(980) 


O . \) 1 -J— \J . i *J 


1 64 


0.28 


-0.19 


0.69 


Q1 


0/o(98O) (°) 


156±3 


41 


-2.2 


4.3 


-0.78 


29 


0^.(1020) (a.u.) 


1.15±0.02 


-0.02 


0.04 


0.003 


0.06 


-0.01 


^0(1020) (°) 


-15±4 


32 


-13 


0.6 


-10.4 


26 


a /o(1710) ( a - u -) 


1.27±0.07 


-0.83 


0.06 


-0.07 


0.22 


-0.87 


^/o(1710) (°) 


102±4 


-27 


-9.4 


3.0 


-6.7 


-15 


Oadd (a.u.) 




5.2±0.4 


1.77±0.21 0.92±0.15 


6.3±0.9 


2.27±0.17 


4>a,dd (°) 




193±4 


93±6 


-179±16 


117±9 


51±4 


X 2 /v 


278/119 


192/117 


249/117 


241/117 


256/117 


200/117 



TABLE VI: Fits to CLEO-c data using the E687 model with additional K + K resonances. For 
the contributions that do not change the entries in the table are changes from the E687 model. 



Parameter 


E687 Model 


/ 2 (1270) 


a 2 (1320) 


/o(1370) 


/o(1500) 


/ 2 (1525) 


a (1450) 


0(1680) 


m K*(892) 


895.8±0.5 


-0.4 


-0.1 


-0.9 


-0.5 


0.0 


-0.8 


0.1 


rV*(892) 


44.2±1.0 


2.3 


2.4 


1.5 


0.6 


0.6 


1.0 


1.2 


a ^(1430) ( a - u 


1.76±0.12 


0.11 


0.08 


-0.25 


-0.03 


-0.16 


-0.22 


-0.18 


0K*(143O) (°) 


145±8 


-32 


-28 


1.0 


-15 


1.7 


-15 


18 


a/o(980) ( a -u-) 


3.67±0.13 


0.29 


0.26 


1.05 


0.52 


0.03 


1.09 


0.20 


<^/ (980) (°) 


156±3 


-2 


-1.6 


1.3 


2.3 


0.22 


3.8 


10.5 


a 0(lO2O) ( a - u 


1.15±0.02 


-0.03 


-0.04 


-0.02 


-0.003 


-0.02 


-0.007 


-0.012 


<?V(1020) (°) 


-15±4 


-7 


-6.3 


7.2 


-0.6 


1.5 


4.3 


13.2 


«/o(1710) ( a - u -) 


1.27±0.07 


0.08 


0.07 


-0.16 


0.17 


-0.04 


0.03 


-0.018 


0/o(171O) (°) 


102±4 


7 


4.7 


-13 


-4.1 


-3.8 


-17 


5.3 


Oadd ( a -U.) 




0.64±0.09 0.45±0.06 1.15±0.09 0.50±0.05 0.50±0.07 1.32±0.10 1.04±0.17 


</>add (°) 




17±9 


40±8 


53±5 


132±7 


173±10 


103±5 


-4±11 


X 2 /" 


278/119 


237/117 


237/117 


178/117 


229/117 


249/117 


192/117 


256/117 



events for data plus the signal MC efficiency or background box to constrain the variation 
of the efficiency or background parameters. We also fit allowing the signal fraction to float, 
and find / S j g = 0.8495 ± 0.0070 which is consistent with 0.8490 used in the central fit. 

We estimate a systematic uncertainty of the Model A fit parameters by combining the fit 
results from Tables IVHt IVIIIt |X| Kit and Kill None of the systematic variations dominate 
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1 1.5 2 2.5 3 0.5 1 1.5 2 



m 2 (K + K ) (GeV/c 2 ) 2 m 2 (K-7i + ) (GeV/c 2 ) 2 

3541 208-023 

3541208-025 I 




m 2 (K + 7t + ) (GeV/c 2 ) 2 m 2 (K + K") (GeV/c 2 ) 2 



FIG. 4: Fit to data for Model A, and projections of the Dalitz plot. The final plot shows the 
m 2 (KK) projection of Dalitz plot for values of m 2 (KK) larger than the contribution from the 
0(1020). 

the uncertainty. The systematic uncertainty is estimated as the mean change from the 
central fit result, 5 Mean, added in quadrature to the RMS of all variations. The resulting 
systematic uncertainties on the parameters are given in Table IXIIII 

V. CONCLUSION 

We perform a Dalitz plot analysis of the — > K + K~it + decay with the CLEO-c data 
set of 586 pb _1 of e + e _ collisions accumulated at yfs = 4.17 GeV. This corresponds to about 
0.57 million DfD*~ pairs from which we select 14400 candidate events with a background 
of 15%. We compare our results with the previous measurement from E687 using the isobar 
model and find good agreement with the E687 parameters, as shown in Table IIVI We find 
that all resonances from E687 model are significant and their exclusion degrades the fit 
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TABLE VII: Fits to data using Model A with additional non-resonant or K + tt resonance. For 
the contributions that do not change the entries in the table are changes from Model A. 



Pqvq T~yi 0"f" or 
JT dl cLLUC LCI 


Mnrlpl A 

1VXUU.C1 -TV 


MR 

iY it 


if* (1410) #1(1430) if* (1680) 


K' 
r\/ 


m K*(892) 


OQ/I Q-|_n c; 

oy^t.yznu.o 


U.O 


n i 

U.l 


n 9 
u.z 


n 9 
u.z 


n 1 

-U. 1 


1 X*(892) 


7-1-1 1 


n i 

-U. 1 


n i 

U.l 


n 8 
U.o 


U.D 


-U.O 


a i^(1430) (a.u. J 


i _i_n 1 1 


n 1 878 
-U. lo / o 




U.UUUO 


n i a^a 

-U. L^kO^k 


U.ZDoO 


9K*(U30) ( ) 




i n 8^^ 


P. A A Ad 

U.4440 


-4.8755 


2.3676 


- i .OOUo 


a /o(980) ( a - u 


4. / z±U.lo 


-U.UOzy 


0.0U57 


0.2566 


-0.2530 


-0.2078 


^/o(980) ( J 


lb / ±o 


o ii cro 


-0.7457 


1.0875 


1.7545 


-4.5506 


00(1020) ( a - u 


i i Q_i_n no 
l.lo±U.0z 


-U.0005 


-0.0001 


-0.0096 


-0.0159 


0.0047 


^(1020) I ) 


8-1-4 


"i QQ7^ 
o.yy / o 


-0.1144 


-4.8349 


5.2172 


-5.0235 


a /o(1370) V a - U -J 


i i cr_i_n no 
i.ioznu.uy 


-u.uy / y 


-0.0055 


0.0535 


0.0103 


0.0890 


<?>/o(1370) I J 




U.UtJUU 


-1.6829 


-4.4427 


3.2688 


-11.386 


a /o(1710) V a - U d 


iii _i_n n7 

l.llznU.U i 


-U. 10UZ 


-0.0093 


-0.0157 


-0.0442 


-0.0940 


?Mmo) [ ) 


804-^ 
oy in (j 


7 ^1 9fi 


-1.2087 


3.7678 


2.4526 


-6.2195 


«add (a.U.) 


u 


i.oztU.D 


0.10±0.13 


1.00±0.26 2.18±1.33 0.50±0.18 


<Pa,dd { ) 


n 
u 


1 /1 7-1-1 Q 


-3±119 


105±11 


-72±13 


163±25 




n 

U 


1 ^4-1 4 


0.01±0.03 0.40±0.22 0.30±0.44 0.40±0.32 


FF add (%) @ 90% C.L. 





<3.3% 


<0.05% 


<0.7% 


<0.9% 


<0.8% 


FF[if*(892)] (%) 


47.4±1.5 


47.5 


47.5 


47.8 


48.3 


47.5 


FF[if *(1430)] (%) 


3.9±0.5 


3.0 


3.8 


4.4 


3.3 


5.5 


FF[/ (980)] (%) 


28.2±1.9 


27.7 


28.4 


32.3 


26.2 


25.7 


FF[0(1O2O)] (%) 


42.2±1.6 


41.9 


42.1 


42.3 


42.1 


42.1 


FF[/ (1370)] (%) 


4.3±0.6 


3.5 


4.2 


4.8 


4.5 


4.9 


FF[/ (1710)] (%) 


3.4±0.5 


2.6 


3.4 


3.4 


3.3 


2.9 


FF R (%) 


129.5 


127.8 


129.4 


135.4 


127.9 


129.0 




178/117 


174/115 


177/115 


170/115 


175/115 


173/115 



quality. 

However, the fit quality is signficantly improved if we add an additional K + K~ reso- 
nance to the model. As shown in Tables |V] and IVIl almost any additional resonance or 
non-resonant contribution improves the agreement with the data. The best improvement is 
achieved if we add an / (1370)tt + contribution. We find that a six-resonance model, con- 
taining contributions from K*(892) K + , K* (1430) K + / (980)7r+, 0(1O2O>+, / (1370)tt+ 
and /o(1710)7r + resonances, gives better consistency with our data with x 2 jv = 178/117. 
Tables IVHI and IVIIII show that any further additional resonance does not have a significant 
amplitude, fit fraction, or significantly improve the fit quality and we give upper limits on 
their fit fractions at the 90% C.L. 

In Table HXl we show the resonance parameters when they are allowed to float in the fit. 
We find that the if*(892) width is 5 MeV/c 2 smaller than in PDG. This result is consistent 
with our observation in the D + — > K~ir + ir + analysis [l5|. Other resonance parameters are 



consistent with their values from the PDG [8] or the BES experiment [13| for / (980) 



We estimate a systematic uncertainty on fit parameters from numerous fit variations, 
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TABLE VIII: Fits using Model A with additional K + K resonance. For the contributions that 
do not change the entries in the table are changes from Model A. 



Parameter 


Model A 


/ 2 (1270) 


a 2 (1320) 


/o(1500) 


/ 2 (1525) 


a (1450) 


0(1680) 


m K*(892) 


894.9±0.5 


-0.5 


-0.3 


-0.1 


-0.2 


-0.1 


-0.1 


r_ft"*(892) 


45.7±1.1 


1.2 


1.2 


-0.1 


0.2 


-0.1 


-0.2 


a ftT*(1430) ( a -U.) 


1.51±0.11 


-0.0518 


-0.0587 


-0.0060 


-0.0822 


-0.0210 


-0.0152 


01^(1430) (°) 


146±8 


-13.610 


-7.5258 


1.1483 


0.1662 


2.4740 


-0.8833 


a fnfQ80 , l (a.u.) 


4.72±0.18 


0.0864 


-0.0037 


0.0521 


-0.0239 


0.1123 


0.0113 


(yfn Cqsri'i (°) 


157±3 


-0.6746 


-0.6856 


0.6617 


-0.3009 


1.1151 


-0.1360 


°</>(i020) (a.u.) 


1.13±0.02 


-0.0105 


-0.0126 


0.0058 


-0.0058 


0.0068 


0.0056 


00(1020) (°) 


~8±4 


-2.1292 


-1.5385 


0.5046 


-0.1244 


1.2202 


-0.4788 


«/o(1370) (a-u.) 


1.15±0.09 


-0.0176 


-0.0343 


0.0336 


-0.0168 


0.0150 


-0.0039 


0fn(137O) (°) 

J U \ J. 


53±5 


1.0892 


-0.3964 


3.8125 


1.4021 


14.6004 


0.3390 


a fn(1710) (a.U.) 


1.11±0.07 


0.0041 


-0.0165 


-0.0161 


-0.0100 


-0.0533 


0.0007 


0fnd71O) (°) 


89±5 


4.7785 


2.7846 


-1.9584 


-2.2626 


-3.6665 


-0.9276 


Oadd (a.U.) 





0.40±0.09 0.26±0.06 0.07±0.04 0.23±0.08 0.37±0.28 


0.10±0.16 


0add (°) 





22±14 


51±15 


37±66 


180±26 


24±17 


-93±122 


FF add (%) 





0.24±0.11 0.20±0.09 0.04±0.10 0.09±0.05 0.38±0.60 0.008±0.031 


FF add (%) @ 90% C.L. 





<0.4% 


<0.3% 


<0.17% 


<0.16% 


<1.2% 


<0.05% 


FF[K*(892)] (%) 


47.4±1.5 


47.2 


47.4 


47.3 


48.0 


47.3 


47.4 


FF[K *(1430)] (%) 


3.9±0.5 


3.8 


3.7 


3.9 


3.6 


3.8 


3.8 


FF[/ (980)] (%) 


28.2±1.9 


30.0 


29.0 


28.8 


28.4 


29.4 


28.2 


FF[0(1O2O)] (%) 


42.2±1.6 


42.1 


42.2 


42.2 


42.1 


42.2 


42.1 


FF[/ (1370)] (%) 


4.3±0.6 


4.2 


4.1 


4.5 


4.2 


4.3 


4.2 


FF[/ (1710)] (%) 


3.4±0.5 


3.5 


3.4 


3.3 


3.4 


3.1 


3.4 


Eh FF R (%) 


129.5 


131.1 


130.2 


130.0 


129.8 


130.5 


129.3 


x 2 /^ 


178/117 


169/115 


170/115 


177/115 


172/115 


176/115 


178/115 



and Table [XTTTI shows the final results on fit parameters with their statistical and systematic 
uncertainties. 
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TABLE IX: Optimal resonance parameters. The uncertainties for the CLEO-c results are statistical 
only. 



Resonance Parameter (MeV/c 2 ) 


Ppntrnl Fit 


Flnntprl 


pnr r»i 

L 1 JVl IOI 


Pi. [OilZ ) 


in 


cnc; S-l-n ^ 
oyo.oztu.o 


one; S-i_n ^ 


qqr nn-i-n 9^ 
oyo.uuztu.zo 




p 
i 


/m 9_i_i n 


A A 9_|_1 n 


Kn Q-i_n R 


A Q (14,111; 


7?? 


1/11/1 


1 /199_|_9Q 
14ZZztZo 


1/11 /i-i-f; 

1414k ItO 




r 


9Q0 
zyu 


zoy zn^±o 


9Qn.4-91 
zyuznz i 


/o(980) 


m 


965 


933±21 


980±10 




QlTTT 


406 


393±36 


r=4o to ioo 




9KK 


800 


557±88 




0(1020) 


m 


1019.460 


1019.64±0.05 


1019.460±0.019 




r 


4.26 


4.780±0.14 


4.26±0.05 


/o(1370) 


m 


1350 


1315±34 


1200 to 1500 




r 


265 


276±39 


200 to 500 


/o(1710) 


m 


1718 


1749±12 


1718±6 




r 


137 


175±29 


137±8 
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TABLE X: Fits to data using Model A with floating resonance parameters. After the first column 
of data the entries in the table are changes from Model A when the parameters of resonance at the 
top of the column are allowed to float. 



Parameter Model A K*(U30) / (980) 0(1020) / (1370) /o(1710) 



m K*(892) 


Q04 q_i_n c 


-D 1 


n 


n 9 


-D 1 


n 1 


A (892) 


45.7±1.1 


-0.1 


0.2 


0.1 


0.0 


-0.5 


ft, r^*/-iyion\ (311 1 


1.51±0.11 


-0.1449 


-0.1527 


0256 


0533 


-0 0305 


< PA"q(1430) I ) 


1 4«+Q 

J.'iU—l—O 


8 6060 


-3 2558 


10 2102 


7 5225 


-5 6685 


a /o(980) (,a.U.) 




U.UtJ / u 


u.oo / 


u.ou / 


n 0^40 

U.UO^iU 


1 7R7 

u. ± 1 u / 


^/o(980) (°) 




1 1 909 


1 ^ £84 


n 00^7 


1 9907 


^ 40£8 


"0(1020) l^d.U.J 




U.UUJO 


n nni s 


U.U ( ou 


n 00^7 

U.UUO 1 


n m R7 


^(1020) (°) 


-8+4 


-D 891 6 

VJ . O Zv ± u 


5 99Q1 


1 56Q7 


n Q61 3 


1 3374 


u /o(1370) 


1 1 K_i_n no 

± . ± <J_I_U . \JU 


D 0473 


-0 031 Q 


-D 0508 


D D9Q3 


-0 1948 


</>f n (4370 , | (°) 


53±5 


-2.5387 


4.8538 


-2.6304 


-17.247 


3.0673 


a /o(1710) ( a - u -) 


1.11±0.07 


-0.0060 


-0.0096 


-0.0291 


-0.0656 


0.4223 


0/o(171O) (°) 


89±5 


-1.9306 


-1.2058 


-2.4148 


0.0913 


20.0144 


FF[K*(892)] (%) 


47.4±1.5 


47.3 


47.2 


47.4 


47.5 


46.8 


FF[K *(1430)] (%) 


3.9±0.5 


3.8 


3.2 


4.1 


4.2 


3.7 


FF[/ (980)] (%) 


28.2±1.9 


27.5 


29.7 


24.8 


27.7 


29.7 


FF[0(1O2O)] (%) 


42.2±1.6 


42.2 


41.8 


43.3 


42.2 


42.0 


FF[/ (1370)] (%) 


4.3±0.6 


4.6 


4.0 


3.9 


4.4 


3.3 


FF[/ (1710)] (%) 


3.4±0.5 


3.4 


3.4 


3.3 


3.0 


4.1 


Er ff r (%) 


129.5 


128.8 


129.2 


126.8 


129.0 


129.5 




178/117 


177/115 


169/114 168/115 176/115 166/115 
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TABLE XI: Fits to a variety of data samples using Model A with central efficiency and background. 
After the first column of data the entries in the table are changes from Model A with the variation 
indicated at the top of the column. 



Variation 


Central 


Fit 


Early 


Late 


Only 


Only 


Tight Loose 


Low Side 


High Side 


Parameter 


Model 


A 


Data 


Data 


D+ 


DJ 


la X la 3a X 3a 


Band 


Band 



m if*(S92) 
IV* (892) 
a if* (1430) ( a - u -) 
0if*(143O) (°) 
a /o(980) ( a - u -) 
0/ o (98O) (°) 
a ^(l()20) (a.u.) 
00(1020) (°) 
O/ (1370) ( a - u -) 
0/0 (1370) (°) 
a /o(1710) ( a - u -) 



894.9±0.5 
45.7±1.1 
1.51±0.11 

146±8 
4.72±0.18 

157±3 
1.13±0.02 

-8±4 
1.15±0.09 

53±5 
1.11±0.07 

89±5 



-0.4 

0.1 
0.0138 
-10.971 
0.3277 
-1.3604 
0.0053 
-2.9134 
0.0976 
-2.8318 
0.0786 
-3.3881 



3.0 

0.0 
0.0177 
9.7985 
-0.3513 
1.1808 
0.0008 
2.2119 
-0.1031 
2.2204 
-0.0830 
2.2247 



-0.7 
-0.8 
-0.0023 
-17.161 
0.0484 
-6.3697 
0.0084 
-8.3156 
-0.0131 
-4.6088 
-0.0412 
0.1313 



0.7 

0.8 
0.0398 
17.257 
-0.0416 
6.6295 
0.0011 
8.5410 
0.0250 
2.4167 
0.0483 
-0.5966 



-0.2 
-0.2 
-0.0205 
-6.2148 
-0.0364 
-4.5506 
0.0153 
-7.0696 
-0.1193 
-6.5716 
0.0403 
0.7797 



0.2 

1.0 
-0.1276 
14.408 
0.0244 
2.8515 
-0.0049 
5.5073 
0.1395 
-1.3470 
0.0070 
2.6467 



-1.2 

4.8 
-0.8084 
18.400 
0.6752 
3.1875 
0.0079 
8.5766 
0.5111 
-14.394 
0.1877 
16.146 



-1.4 

2.2 
0.7309 
-66.057 
0.3610 
-23.699 
0.0210 
-35.140 
0.3938 
-28.267 
-0.3847 
-5.0150 



FF[i<"*(892)] (%) 


47.4±1.5 


47.2 


47.7 


47.9 


46.7 


47.2 


46.8 


43.4 


48.0 


FF[E"*(1430)] (%) 


3.9±0.5 


4.0 


4.1 


3.9 


4.2 


3.8 


3.3 


0.9 


8.3 


FF[/ (980)] (%) 


28.2±1.9 


32.1 


24.4 


28.6 


27.9 


27.6 


28.8 


37.6 


31.5 


FF[</>(1020)] (%) 


42.2±1.6 


42.0 


42.3 


42.1 


42.2 


42.7 


42.0 


43.3 


41.8 


FF[/ (1370)] (%) 


4.3±0.6 


5.0 


3.5 


4.1 


4.5 


3.4 


5.4 


9.0 


7.4 


FF[/ (1710)] (%) 


3.4±0.5 


3.9 


3.0 


3.2 


3.8 


3.7 


3.5 


4.8 


1.4 




129.5 


134.2 


124.9 


129.7 


129.2 


128.3 


129.9 


138.9 


138.4 




178/117 


134/117 203/117 166/117 123/117 155/117 201/117 


140/117 


138/117 


Events on DP 


14400 


7334 


7066 


7233 


7167 


7200 


19177 


6682 


7232 


/sig 


0.8490 


0.8518 


0.8466 


0.8496 


0.8497 


0.9238 


0.7484 


0.4338 


0.5696 



19 



TABLE XII: Fits to data using Model A with floating efficiency and background coefficients, fits 
with floating / s i g , and with floating background coefficients Bk* and for the narrow resonance 
contributions to the background. After the first column of data the entries in the table are changes 
from Model A with the variation indicated at the top of the column. 



Parameter 


Model A 


Float Ei 


Float Bi 


Float U\cr 


ln K*(S92) 


894 q+n 5 


o 


0.1 


o 




45.7+1.1 





-0.2 


o 


CLis* fi /iqr^ ( EL U. ) 


1.51+0.11 


-0.0018 


-0.0121 


0.0023 


(\~) 7>" * / i ion'\ ( I 

Srif J (1430 J V / 


146+8 


0.1630 


-1.6971 


0.2116 


a / (980) ^a.u.j 




-0 0096 


-0 0339 


-0 0043 


^/o(980) I ) 


157+3 


3362 


-0 6851 


2704 


a ^(1020) V* 1 - u * ) 


1 13+0 D2 


0.0034 


-0.0007 


0028 

\J • \J\J jLJKJ 


c r</)(1020) I / 


-8+4 


0.1282 


-0.9907 


-0.0391 


u /o(1370) V a - U -^ 


1 15+0 09 


-0.0015 


0.0112 


0006 

\J • \JKJ\JKJ 


0fnd37(u (°) 


53+5 


0.1323 


-0.5403 


0.0792 


a /o(1710) ( a - u 


1.11±0.07 


-0.0007 


-0.0539 


-0.0038 


0/o(171O) (°) 


89+5 


-0.2072 


-1.1088 


-0.3882 


FF[if*(892)] (%) 


47.4+.1.5 


47.4 


47.7 


47.4 


FF[tfg(1430)] (%) 


3.9+0.5 


3.9 


3.9 


3.9 


FF[/ (980)] (%) 


28.2+1.9 


28.2 


28.1 


28.2 


FF[^(1020)] (%) 


42.2+1.6 


42.2 


42.2 


42.2 


FF[/ (1370)] (%) 


4.3+0.6 


4.2 


4.4 


4.3 


FF[/ (1710)] (%) 


3.4+0.5 


3.4 


3.1 


3.4 


E« FF^ (%) 


129.5 


129.4 


129.3 


129.4 


X 2 /v 


178/117 


679/562 270/188 


178/116 
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TABLE XIII: Summary of systematic cross checks for Model A. Fit parameters are shown with 
their statistical and systematic uncertainty respectively. The "JMean" and "RMS" account for 
variation of the fit parameters in the systematic cross checks as discussed in the text. The "Total" 
is a quadratic sum of "JMean" and "RMS" and after rounding is the systematic uncertainty given 
in the second column. The results of the E687 Model are also shown for comparison. 



P stq m PT' 


]\/r nr ip1 A 


JMean RMS Total E687 Model 


m K*(&92) 


004 q+n 5+n 7 


0.088 


0.654 


0.660 


895 8±0 5 


1 i<T*(892) 


4c 7_i_i 1 _i_n c 


0.148 


0.499 


0.520 


44.2±1.0 


U _K*(1430) V*- U V 


1 tci _i_n 1 1 _i_n nn 

X • O J l_ Ui 1J 1— \J • \J *J 


-0.024 


0.089 


0.092 


1 76±0 12 


<Pi<"*(1430) I ) 




-0 693 


8 449 


8 465 


1 4.c:_i_q 


a /o(980) ( a - u 


A 794-fl 1 84-fl 1 7 




U. 1U 1 


n 1 70 


o.u / zcu. 10 


0/ o (98o) I J 


1 c:7_i_o_i_ / i 




4 n^fi 

'i.UOU 


A DM 




°0(iO2O) (a.u.) 


1 1 Q-i_n n.9-i-n no 

l.lOItU.UZltU.UZ 


Odd 


01 7 


yj. \j ±0 


1 1 c: 1 n no 

JL . JLfJ -i—\J ,\J 


<?V>(1020) I J 


_0_|_/|_|_/| 


081 


3 850 


3 851 

*J. (J) (J X 


-15±4 


«/o(1370) l^a.u.j 


1 1 K_i_n fiQ+f) Dfi 
± . ±(jznu.uyznu.uu 


-0.003 


0.063 


0.063 




^/o(1370) (°) 


53±5±6 


-0.536 


5.820 


5.845 




a /o(mo) (a-u.) 


1.11±0.07±0.10 


-0.004 


0.098 


0.098 


1.27±0.07 


0/o(171O) (°) 


89±5±5 


0.195 


4.916 


4.920 


102±4 


FF[if*(892)] (%) 


47.4±1.5±0.4 


0.016 


0.357 


0.4 


48.2±1.2 


FF[i^(1430)] (%) 


3.9±0.5±0.5 


0.036 


0.460 


0.5 


5.3±0.7 


FF[/ (980)] (%) 


28.2±1.9±1.8 


0.096 


1.792 


1.8 


16.8±1.1 


FF[0(1O2O)] (%) 


42.2±1.6±0.3 


0.018 


0.277 


0.3 


42.7+1.3 


FF[/ (1370)] (%) 


4.3±0.6±0.5 


0.044 


0.488 


0.5 




FF[/ (1710)] (%) 


3.4±0.5±0.3 


0.044 


0.311 


0.3 


4.4±0.4 


Eij FF R (%) 


129.5±4.4±2.0 


0.020 


1.981 


2.0 


117.3±2.2 


x 2 /^ 


178/117 








278/119 
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